Reliability of the Heitler-London approach for the exchange coupling 
between electrons in semiconductor nanostructures 
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We calculate the exchange coupling J between electrons in a double-well potential in a two- 
dimensional semiconductor environment within the Heitler-London (HL) approach. Two functional 
forms are considered for the single-well potential. We show that by choosing an appropriate and 
relatively simple single-electron variational wave function it is possible, within the HL approach, to 
significantly improve the estimates for J. In all cases the present scheme overcomes the artifacts and 
limitations at short interdot distances, previously attributed to the HL method, where unphysical 
triplet ground states have been found, and leads to an overall agreement with analytic interpolated 
expressions for J obtained for a donor-type model potential. 

PACS numbers: 03.67.Lx, 85.30.-z, 73.20.Hb, 85.35. Gv, 
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For quantum computer architectures based on elec- 



tron spins, 



1,2,3,4,5 



two-qubit logical gates can be imple- 



mented by entangling electrons bound to neighboring po- 
tential wells through a transient exchange coupling^ Ac- 
curate evaluation of this coupling — namely, the differ- 
ence in energy between triplet and singlet two-electron 
ground-state configurations — requires numerically inten- 
sive calculations which are highly dependent on the phys- 
ical parameters of the system. Such parameters are in 
general not precisely known, and for this reason the 
Heitler-London (HL) method^ 7 - is of special interest in 
providing intuitive insight and computational simplic- 
ity for prospective estimations of the exchange coupling 
strength over a wide range of potential-related param- 
eters. This is particularly important for gated quan- 
tum dots in which the exact form of the potentials, pro- 
duced by electrodes over a two-dimensional electron gas 
(2DEG), has to be modeled. Also, within the HL ap- 
proach, it is in principle possible to explore the exchange 
energy at large interwell distances, which would require 
great computational effort in more accurate methods. 

Prospective model calculations based on the HL 
method have led to an inaccurate exchange coupling 
asymptotic (large interdot distances) decay and to un- 
physical artifacts such as predicting a triplet ground state 
at short interwell distances (negative- J anomaly). 8,9 
It is usually argued that these problems are due to 
limitations in the HL approach^* 9 - and it has been 
suggested that they could only be overcome within 
more complex quantum chemistry or numerical meth- 
ods ,12,1^12,13,14,15 .mm* We show here that the limita- 
tions of previous studies are mainly due to the choice 
of the single-particle wave functions rather than being 
intrinsic to the HL method. 

Within the single-valley effective mass approximation, 
the Hamiltonian for an electron trapped in a single cylin- 
drically symmetric well potential in a 2D semiconduc- 
tor environment is given by h(p) — + V(p), where 



P = \J x2 + y 2 is the distance from the well minimum, 
and the energies and distances are given in atomic units 
Ry* = m^e 4 /2h 2 e 2 and a* — H 2 e/m^e 2 , respectively. 
Here the semiconductor is characterized by m±, the 
transverse effective mass, and by the dielectric constant 
e. For Si, the experimental values of m± and e lead to 
a* Si = 3.157 nm and Ry^ = 19.98 meV. For GaAs-based 
quantum dots, these values are a* GaAs = 10.343 nm and 



As 



5.31 meV. The first term in the Hamiltonian 
h refers to the kinetic energy of the electron, and the sec- 
ond term is the model potential, which depends on the 
architecture of the physical system involved. 

The confining potential produced by a simple shallow 
donor (e.g., P) positioned a distance d from an interface 
with a barrier is given by 



V(p) 



(i) 



This does not include the effect of the image charge 
in the Si/Si02 architecture, which could be accounted 
for straightforwardly by multiplying the potential by 
2esi/(esi + £Si0 2 ) ~ 1-5. Near the potential minimum 
we get V(p » 0) = —2/d + p 2 /d 3 , defining the parabolic 
approximation for this model potential. 8 

For gated quantum dots, several model potentials have 
been adopted in the literature. In particular, the func- 
tional form in Eq. ([T]) yields a nearly harmonic confine- 
ment around the dot center and asymptotically follows 
a Coulombic decay at large distances from the center 
[V(p — ► oo) rs —2/p]. Previous studies of gated quan- 
tum dots (e.g., Refs. [ililfllllTil and fl9h propose a faster 
decay for the quantum dot confining potential, given by 
a Gaussian form 



V'(p) = ~V e-P 2 / L2 



(2) 



We also consider this model potential below. 

For d = the Hamiltonian with V(p) in Eq. ([1]) reduces 
to that of a 2D hydrogen atomj&S . so the exact ground- 
state solution is d> 



evo exp(— aop)/v27r, with a, 
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FIG. 1: (Color online) Ground-state energies obtained for 
the three approximations considered here. As the dot-dashed 
lines, the much larger energies obtained through the parabolic 
Gaussian (PG) scheme. As the dashed lines, the energies 
for the variational Gaussian (VG), which are very close to 
those calculated through the more general matched varia- 
tional (MV), given by the crosses. The insets show the percent 
reduction on the energy from VG to MV. Results in (a) corre- 
spond to the potential V(p) in Eq. Q while results for V'(p) 
in Eq. ([2]) are given in (b). 



0.5a*, and the energy is E(d = 0) = -4Ry*. In the gen- 
eral case (d ^ 0), some approximation must be adopted. 
The simplest approach consists of using the exact ground- 
state wave function of the parabolic potentials defined in 
the paragraph following Eq. |T]). We call this approach 
the parabolic Gaussian (PG), and the ground state has 
the functional form (f)pc — ftpa/V^ ex P (~PpgP 2 /?) , 
where (3pc — (2/d 3 ) 1 / 4 , and the energies are given by 
Epq = 2(3 2 — 2/d. A better approach^! is to treat (3 as 
a variational parameter [variational Gaussian(VG) wave 
function approach] keeping the exact form of the poten- 
tial. This leads to significant energy reduction as com- 
pared to PGj2^ as shown in Fig. DJa). 

The VG form of the wave function is well suited close 
to the potential minimum (where the parabolic fit is more 
appropriate) . At large distances from the well minimum, 
however, the potential is finite and we expect the ground- 
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FIG. 2: (Color online) (a) Illustration of \4>mv\ for d — 6a* 
[see Eq. J3)]. (b) The optimized variational parameters for 
MV, fj, and 1/(3, are shown here to be of the same order of 
magnitude. The given results correspond to a single well of 
the form V(p) as in Eq. (p}. 



state wave function to decay exponentially rather than as 
a Gaussian. The exponential decay is expected for any 
electronic bound state (Eq) in a potential that goes to 
zero at large distances, contrary to the Gaussian decay, 
characteristic of the infinitely confining parabolic well. 
This argument also applies to the Gaussian model po- 
tential V'(p) in Eq. @. To correctly take into account 
the form of the potential at short p — > and long p — > oo 
distances, we define a variational wave function by parts: 
(i) a Gaussian behavior up to some distance p < /i, suit- 
able for a nearly parabolic attractive potential, and (ii) 
an exponential decay for p > /i, adequate for a finite po- 
tential. This matched variational (MV) wave function 
[illustrated in Fig. [2ja)] is written 



(P) 



fA ie xp(-^) if 
\^ 2 exp(-^) if 



P < /!, 
p > (J,, 



(3) 



and involves the parameters A±, A2, (3, a, and /1. From 
the continuity of the wave function and of its first deriva- 
tive at p = fj,, and normalization, two independent vari- 
ational parameters are left, which we choose to be (3 and 
p,. These are obtained (Fig. ^ by minimizing the energy 
expectation value. The parameter a is obtained through 
the constraint a = (3 2 p. Note that the MV wave function 
reproduces VG for 1/(3. 

The variational principle regards only the value of the 
energy, drawing no conclusions on the validity of the trial 
wave function.— Since the main contribution to the ex- 
pectation value of the energy is given by the p values 
closer to the potential minimum, the more general MV 
does not lower the energy much (typically less than 1%) 
with respect to the VG, as shown in the insets in Fig. [T] 

In general, the optimized variational parameters p, and 
1/(3 are of the same order of magnitude, as shown, for 



3 



example, in Fig. [2jb). This means that, although the 
energy gain of MV with respect to the VG is not very 
significative, the form of the wave functions is different, 
particularly at the tails. The functional form of the wave 
function decay is relevant for the exchange coupling, as J 
depends very strongly on the neighboring wave-function 
overlap. Around the potential minimum, however, the 
difference between the VG and MV wave functions is not 
so important, since Pvg ~ Pmv f° r d > 3a*. 

We write the Hamiltonian for a two-electron system in 
a double well as 

H(l,2)= [-V 2 l +V DW (z)}+—, (4) 

where i labels each electron and the last term is the 
electron-electron Coulomb repulsion. The double-well 
potential Vhw(i) consists of identical wells a distance 
R from each other, 

V D w{i) = V (xi - ^,y^j +V \Xi + ^Vij ■ (5) 

We denote as \cf>A{i)) and |</>s(z)) the one-electron or- 
bitals centered at +R/2 and —R/2, respectively. Within 
the HL approach, the exchange coupling is calculated as 
J = Etripiet — E S i„gi et with the singlet and triplet wave 
functions written as the symmetric and antisymmetric 
combinations of the one-electron orbitals^ 

Anomalies in the HL approach within the PG and VG 
have been pointed out in the literature. In the presence of 
a finite interwell barrier, J should decay exponentially 24 
for R — > oo. This cannot be obtained through PG or VG 
orbitals, as both have a Gaussian decay at all distances. 
Besides the failure in reproducing the expected asymp- 
totic behavior, the PG and VG also lead to unphysical 25 
negative values of the coupling at small R&2- 

In Fig. [3l we show the calculated exchange coupling 
for the three forms considered for <j>. We also show the 
trend of the dependence of J upon the geometrical pa- 
rameter d by showing the results for d = 2a*, 4a*, and 
6a*. There are two main results to stress here: (i) while 
the PG and VG approaches give unphysical negative val- 
ues for the exchange at short distances R, the MV wave 
function gives a positive J for the entire range of param- 
eters studied^ (ii) at large i?, MV gives the expected 
exponential decay for J(R) (as inferred from the linear 
behavior in the log scale of the figure). This decay is 
much slower than the one for VG and PG, leading to 
differences of many orders of magnitude in J among the 
different approaches. 

The "out-of-plane donor" type of potential we are con- 
sidering was previously discussed in Ref. [20|, where an- 
alytic expressions for the exchange integral were inter- 
polated from asymptotic methods which are precise for 
R — ► and R — > oo. This result is compared to ours in 
Fig. [3l where we show that MV is in overall agreement 
with the interpolated expression. 

We expect these results to remain valid for general 
additive double-well models, as long as there is a finite 
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FIG. 3: (Color online) Exchange coupling calculated for the 
indicated values of the geometric parameter d for the quantum 
well potential V(p) in Eq. {T}. The dashed lines and crosses 
are obtained within the HL approach: Dot-dashed lines for 
the PG approximation, dashed lines for VG, and crosses for 
MV. The solid line corresponds to the expression obtained in 
Ref. [13 from the asymptotic behavior for R — > and R — > oo. 
The negative- J anomaly is indicated by the unphysical sharp 
drop in the PG and VG curves at small R values. At large R, 
deviations of PG and VG from the correct asymptotic decay 
increase exponentially. 



barrier between the two wells. Indeed, for the Gaussian 
potential well given in Eq. we find that MV still 
minimizes the energy expectation value [Fig. [fjb)]. The 
calculated exchange coupling, illustrated in Fig.|H shows 
the same trends as obtained for the donor-type potential, 
Eq. |T]), suggesting the general validity of our arguments 
for different types of potential wells. 

The HL approach is a simple tool to estimate ex- 
change coupling between electron spins in artificial and 
real molecules. However, physical inaccuracies reported 
in the past — namely, a triplet ground state predicted for 
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FIG. 4: (Color online) Exchange energies calculated for the 
potential V'(p) in Eq. ©, with V = 0.5Ry* and L = 10a*. 
The general behavior and trends here are equivalent to those 
in Fig. El 



small interdot distances and Gaussian instead of expo- 
nential asymptotic decay at large values of R (Ref. M) — 
severely restrict its applicability. For applications in 
qubits, the small- R behavior of J would be the most 
relevant, since strong enough coupling is needed to guar- 
antee fast gating times. On the other hand, reliable es- 
timates for the large- R coupling must also be available, 
as residual interactions between nominally noninteract- 
ing qubits constitute an additional source of errors affect- 
ing the desired quantum computer operations. Here we 
showed that these shortcomings are not intrinsic to the 
HL approach and may be attributed to the usual Gaus- 



sian form adopted for the one-electron single- well orbital. 
Previous studie a 15 ' 18 establish that, within a Gaussian 
basis set for the single-particle orbitals, overcoming the 
mentioned limitations requires at least 70 molecular or- 
bitals to be considered in solving the two-particle prob- 
lem. We have shown that a variational single-electron 
orbital with appropriate long- and short-distance func- 
tional forms within the HL approach leads to a singlet 
ground state at short R and to the correct exponential 
decay at large R. No quantitative accuracy is attempted 
here, as only the lowest-energy single-electron orbital is 
considered, and even so it is not the exact ground state 
for the single-dot potential. 

Some of us have recently proposed to perform exchange 
operations at the Si/SiC"2 interface between electrons 
laterally bound to donors a distance d from the inter- 
facei 23 ' 27 The form of the potential assumed here ap- 
plies exactly to this problem (except for a multiplicative 
factor related to the image charge, as mentioned previ- 
ously) . The exchange coupling was calculated within the 
HL method through the VG approach for the wave func- 
tion, 2 '' leading to small values of the exchange, < 10 _,i 
meV ~ 10~ 5 Ry*. The results reported here in Fig. [3] 
imply that the actual values of the exchange that can be 
achieved at the Si/Si02 interface are orders of magnitude 
larger and therefore would require much shorter gating 
times than previously expected. 
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